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for  time-harmonic  longitudinal  waves  propagating  in  a one-dimensional 
laminated  composite  with  periodic  structure.  The  development  is  based 
on  the  method  of  Lie-series  and  a hierarchy  of  improved  approximate 
theories  are  obtained  for  wavelengths  longer  than  the  typical  dimensions 


of  the  composite  microstructure.  Effective  material  properties  are 
also  defined  and  the  dispersion  spectrum  compared  with  the  exact  theory, 
over  one  Brillouin  zone.  It  is  also  shown  how  the  effective  proper- 
ties of  the  laminated,  periodic  composite  can  also  be  obtained  by 
homogenization  via  a multiple  scale  expansion. 


Introduction 

When  studying  the  behavior  of  an  inhomogeneous  material  with  material 
properties  which  vary  with  position,  it  often  proves  convenient  to  replace 
it  by  a homogeneous  material,  whose  properties  approximate  that  of  the 
original  material,  in  some  suitable  manner.  The  properties  of  this 
'equivalent'  homogeneous  material  are  then  called  the  effective  or  bulk 
properties  of  the  inhomogeneous  material.  This  idea  is  quite  old  and  can 
be  traced  back  to  Poisson,  Maxwell  and  Rayleigh.  An  extensive  survey 
and  references  on  this  subject  can  be  found  in  the  article  by 
Babuska  [1]*. 

Recently,  with  increasing  use  of  composite  materials  with  laminated 
structure,  and  the  increasing  complexity  of  the  problems  generated  by 
their  study,  the  question  of  developing  'approximate'  theories  using 
effective  properties  has  become  increasingly  important.  Earlier  work  in 
this  area  can  be  classified  in  three  categories,  i)  effective  stiffness 
theory,  ii)  mixture  theory,  and  iii)  low  frequency  long  wavelength 
theory.  Representative  of  these  three  categories  are  the  papers  by 
i)  Sun,  Achenbach,  and  Herrmann  [2],  ii)  Bedford  and  Stern  [3],  and 
iii)  Kohn  [4],  respectively.  In  this  paper  we  outline  in  a compact  and 
formally  exact  form  two  such  theories,  one  using  the  method  of  Lie- 
series  [5]  and  the  other  using  homogenization  via  multiple  scale  expansion. 
We  show  that  the  two  methods  of  developing  the  approximate  theory  for 
elastic  waves  in  laminated  composites  lead  to  equivalent  results.  Further- 
more, in  the  limit  of  low  frequency  and  long  wavelength  these  equations 
are  in  agreement  with  the  equations  obtained  earlier  by  Kohn.  Therefore 
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these  theories,  though  obtained  by  using  three  entirely  different  mathe- 
matical procedures,  may  be  called  long  wavelength  approximation  theories. 
The  accuracy  of  these  approximate  equations  can  be  improved  by  increasing 
the  order  of  che  differential  equation.  Our  results  show  that  the 
approximate  dispersion  spectrum  tends  monotonically  towards  the  exact 
spectrum  by  increasing  the  order  of  approximation  and  the  convergence  is 
from  below.  We  have  shown  that  in  the  case  of  elastic  waves  in  one  dimen- 
sional laminated  composite,  a linear  differential  operator  of  the  form 


O 


l ~ j>  C2 

t -1  11 
n=l 


32n) , reproduces  the  exact  spectrum  in  the  first  Brillouin 


zone  to  a high  degree  of  accuracy. 

The  method  of  Lie-Series  as  used  in  this  paper  has  several  advantages 
over  other  methods  used  in  the  past  in  the  development  of  approximate 
theories.  First  convergence  can  be  justified  and  there  exists  a large 
literature  on  this  subject.  Second,  using  Lie's  exponential  mapping  and 
commutation  theorems  one  gets  the  exact  spectrum  from  the  Lie-form  of 
differential  equations.  The  correspondence  between  a spectral  function  and 
its  differential  equation,  in  conjunction  with  the  implicit  function  theorem, 
then  helps  us  to  develop  a hierarchy  of  improved  approximate  theories. 

Formulation  using  Lie-Series 

We  consider  the  one  dimensional  elastic  wave  equation 

3o/3y  - p32U/3t2  =0,  0<y<“,  t>0  (1) 

o = n3U/3y  , (2) 

where  oeC2  is  the  stress,  UeC2  is  the  longitudinal  displacement,  p the 


mass  density,  n the  stiffness,  y the  longitudinal  distance  parameter 
in  the  undeformed  configuration,  and  t the  parameter  of  time.  The  one- 
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dimensional  medium  is  supposed  to  be  laminated  with  period  2A  = 2(h^+  h^) 
The  unit  cell  is  the  union  of  two  materials;  part  of  the  cell  with  thick- 
ness 2h^  has  material  properties  (p^,n^),  and  the  remainder  of  the 
cell  with  thickness  2h2  has  material  properties  ((>2^2)'  as  shown  in 
Fig.  1.  The  union  of  an  infinity  of  such  unit  cells  gives  us  the  periodic 
laminated  structure  with  period  2A  with  respect  to  y . The  one  dimen- 
sional longitudinal  wave  equation  applies  to  each  part  of  the  cell,  though 
different  parts  of  the  cell  have  different  material  properties. 

The  continuity  of  displacement  and  stress  at  the  interface  A requires 

that 


U^(Y^  + h1,t)  = - h2,t)  . t >_  0 

°1(Y1  + hl,t:)  “ °2(Y2  ' h2,t)’  tiL° 


(3) 


and  similarly  at  the  interface  B,  the  continuity  conditions  are 

Itk/Vk  . , . ITk+l,„k+l  , ^ , . n 

U2(Y2  + h2,t)  *>  (Y1  - h^t) , t >_  0 

k/TJk  , , k+1.  k+1  , . _ 

o2(Y2  + h2,t)  = a1  ( Yx  - h1,t),  tiO 


(4) 


where  the  superscript  k indexes  the  cell.  The  distance  from  the  origin 
of  the  mass  center  of  layer  a(=l  or  2)  of  the  k-th  cell,  1 <.  k < 00 * 

is  denoted  by  Y 

J a 

Equations  (1)  and  (2)  when  combined  can  be  written  as 


(!^  " 6 ft } ^ + 6 ft } U(y»t)  - 0 , 0 y < -,  t > 0 (3) 


where  32  = p/n  is  a constant  and  represents  the  reciprocal  of  the  square 
of  wave  speed.  Since  6 is  a constant,  the  two  differential  operators 
commute  and  the  general  solution  can  be  written  as  the  sum  of  the  solutions 
of  each. 


* 


I- 


u(y,t)  = uo(ct) , t > 0 

y=° 


— W 

where  C = 1/3  and  UQeC  . By  the  exponential  mapping  theorem  the 
solution  of  this  simple  problem  is  given  by  the  Lie-series 

U(y,t)  = eyD  Uo(Ct)  , (7) 

where  the  operator  D is  given  by 

d s e It  • <8> 

Similarly  the  solution  of  the  second  initial  boundary  value  problem 


(12) 
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which  is  the  standard  D'Alembert  solution  of  wave  equation,  and  repre- 
sents two  waves  propagating  with  speed  C without  change  of  form  in 
the  positive  and  negative  directions. 

The  solution  (11)  can  further  be  rewritten  as 

U(y,t)  = (coshyD)UQ(Ct)  + (sinhyD)U1(Ct)  , (13) 

where  = (Uq  + U^)/2  , = (Uq  - U^)/2  and  are  suitable  analytic 

functions  of  the  parameter  t . When  the  operator  D and  the  functions 
UQ  and  are  all  analytic,  there  exists  a positive  number  Y , such 

that  the  Lie-series  (13)  is  uniformly  convergent  for  |y|<  Y , and  there- 
fore termwise  differentiation  is  admissible.  Then  from  Eqn.  (2)  and 
termwise  differentiation  of  Eqn.  (13),  we  find  that  the  stress  is  given 
by 


o(y,t)  = nD[ (sinhyD)Uo(Ct)  + (coshyD) (Ct) ] . (14) 

From  Eqns.  (13)  and  (14)  we  also  find  that 


U(o,t)  - UQ(Ct)  , (15) 

o(o,t)  *=  nDU1(Ct)  * n u(o,t)  h oQ(Ct)  , 

where  li0(Ct)  and  cQ(Ct)  are  both  analytic  functions  of  t , and  repre- 
sent the  values  of  displacement  and  stress  at  y = 0,  respectively.  There- 
fore in  terms  of  Uq  and  oq  the  Lie-series  representation  of  velocity 
and  stress  can  be  written  as 


U(y,t)  = (coshyD)U  (Lt)  + -i-(sinhyD)o  (Ct)  , 
O Up  o 

o(y,t)  = (coshyD) (Ct)  + Cp(sinhyD)Uo(Ct) , 


(16) 


- ^ 


jrm 
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where  dot  indicates  differentiation  with  respect  to  t , and  use  is 
made  of  the  fact  that  the  operator  D commutes  with  the  operators 
coshyD  and  sinhyD  when  and  U^eC111.  We  thus  see  that  if  the  CL 


functions  UQ(Ct)  and  oQ(Ct)  are  the  values  of  displacement  and  stress 


at  some  midpoint  y = 0,  then  in  terms  of  Lie-series  and  exponential 
mapping  theorem,  the  values  of  velocity  and  stress  at  some  typical 
point  (y,t)  are  given  by 


fr(y,t) 

a(y,t) 


coshyD 
Cp  sinhyD 


1_ 

Cp 


vct> 

-°(ct)_ 

(17) 


coshyDj 
0 < y < 

where  coshyD  and  sinhyD  are  Lie  operators  with  the  definitions 


t > 0 


00  00 
C0ShyD  S I M (yD)2n  ’ sinhyD  " I T2n-1)  ! 


(yD) 


2n-l 


(18) 


2 2 

cosh  yD  - sinh  yD  = 1 . 


We  may  remark  in  passing  that  the  matrix  Eqn.  (17)  has  an  inverse,  since 


2 2 

the  determinant  cosh  yD  - sinh  yD  = 1 . 
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Continuity  Conditions  and  Differential  Equation 

At  the  interfaces  A and  B we  must  require  the  continuity  of  stress 
and  longitudinal  velocity.  Thus  from  Eqns.  (3)  and  (4)  the  continuity 
conditions  are 


^1(Y1  + hl’t)  = ^2^2  “ h2,t)  , 


•k,vk  , , *k+l,„k+l  , . 

U2(Y2  + h2,t)  = Ux  (Yl  - h-^t)  , 


t > 0 


°1(Y1  + hi»fc)  = °2(Y2  " h2ft)  ’ 
o2(Y2  + h2,t)  = ok+1(Yk+1  - h1,t)  . 

In  terms  of  midpoint  longitudinal  velocity  (Ct)  and  midpoint  stress 

O f Cl 

oQ  a(Ct)  , <*  = 1,2,  and  the  general  cell  solution  (17),  the  continuity 

conditions  at  the  interfaces  A and  B require  that 


0,1 

+ 

Clpl 

0 i 

0,1 

'S 

Uo,2 

Co, 

2^2 

k 

1 

c 

k 

. k+1 

1 

o,2 

r 

C p 

T 2 

h2 

o ^ 
0,2 

'S 

Uo,l 

C1P1 

, k 

+ 

Co 

s. 

UK 

= C, 

ok  , 

- Cp„ 

0,1 

11 

hl 

0,1 

h2 

o,2 

2 2 

k 

•k 

k+1 

o,2 

+ 

C2P2 

S. 

h2 

U , 
o,2 

= Ch 

hl 

"o,! 

- C1P1 

S.  °n  2 ’ 
h2  0,2 

„ „k+l 
S,  0 , , 


h2  o,2  ’ 


where 


C.  = cosh  hB 
h a a 3t 

a 


S.  = sinh  h 6 — , 
h a a 8t 

a 


B = 1/C  = (p  / n ) 2 , 

a a a a 


a = 1 or  2 


and  a=l,(2)  for  portion  of  the  cell  with  material  properties, 

(^2^2^  anc*  la^er  thickness  2h2,(2h2)  , respectively.  By  first  adding 
and  then  subtracting  the  first  and  subsequently  the  second  pair  of 
equations  we  obtain  the  following  equivalent  set  of  equations 


c,  (Uk+^  + uk  J 

o,l  o,l 

C,  (Uk+|  - Uk  ) 

h^  o , 1 o,l 

sh  ‘"fi  - *1 ,) 

o,l  o,l 

S,  (ik+J  + Uk  ,) 
o , 1 o , 1 


1 . . k+1  k , 

C1P1  hl  °°>1  ' °°.1> 

1 „ , k+1  , k \ 

^lpl  kl  o>l 

1 . . k+1  k . 

^lpl  kl  °’^  °»1 

1 _ , k+1  k . 

Cu  (a  , - o_  .) 


2 c.  IT  . , 

h2  o,2 
2 c „k 

s«  o 


Clpl  hl  0,1  0,1 


C2P2  h2  0,2  ’ 

2 r k 
Clpl  h2  0,2  ’ 

_ 2 s Uk 

ClP]_  h2Uo,2  1 


<22) 


We  shall  now  postulate  the  existence  of  analytic  functions  U^(y,t) 
and  o^Cy.t),  with  respective  radii  of  convergence  greater  than  the  half- 
period A , and  with  the  important  property 


U (Yk  t)  = Uk  (t)  ; 
a a o ,a 


o (Yk  ) 
av  a , t 


k , . 

a . (t) 

o,l 


(23) 


Due  to  the  strong  smoothness  conditions  imposed  on  these  functions  it 
is  possible  to  write 


lim  _A3y, 


Uo,l(t)  = * U!<Y2  " A’C)  = e Vy’C)  * 

Uo!i(t)  = V^1-)  = VY2  + A,t)  = Jiyh  e 3\(y.t)  , 

°o,i^t^  _ °i<‘Yi,t^  = °i^Y2  - A,t^  = y+yj?  6 ai(y,t)  * 

ak^(t)  = o^Y^.t)  = cs^Y^  + A,t)  * y^yk  e ^(y.t)  , 


(24) 


e 
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where  3^  = 9/9y  and  use  is  made  of  Lie's  commutation  theorem.  From 
Eqns.  (24)  and  (22)  it  now  follows  that  at  the  station  the  func- 
tions and  satisfy  the  equations 


li”k  C,  cAu.  - — S,  S o1  - C,  U.  = 0 , 

y^YZL  hi  A 1 CiPi  hi  A 1 h2  2J 

1:L\  C,  S.U..  - — — S,_  C o,  - = 0 , 

y-PY|l_  ^ A 1 C1P1  h1  A 1 C2P2  h2  2_ 

Ch1CA°l  + Ch2°2_  = ° ' 

lim  f • 1 C2P2 

y_>Y2  L hi  A ^ ^l^l  A 1 ^2  2- 

where 


C,  = coshA9„ 
A y 


S , = sinhA9  . . 
A y 


In  addition  we  assume  that  among  the  entire  class  of  C functions 

fitting  the  midpoint  solutions  there  exists  a particular  set  which  sat- 

isfy  Eqns.  (25)  for  all  values  of  y , rather  than  at  the  station  Y2 

only.  This  special  set  of  functions  {U  ,o  } could  then  be  charac- 

a a 

terized  as  a solution  to  the  set  of  four,  infinite  order,  linear  partial 
differential  equations 


C,  CaUi  ~ 7 su  SA°1  = C,  U.  , 

hi  A 1 C]_pi  h^  A 1 h2  2 

C,  M,  ~ 7^—  S.  C o.  = 77^ — S.  c„  , 
h1  A 1 C1p1  h1  A 1 C2p2  h2  2 

S,  Mi  - — C,  C.o,  = - 77^ — C,  0.  , 

h^  A 1 C^p^  A 1 Clpl  ^2  2 


"l^l  “2 


• 1 C2P2 

S.  C.U.  - C,  S a1  = - S.  U,  . 

hj  A 1 Glpl  4 4 1 Clpl  “2  2 
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Now  one  can  express  , and  in  terms  of  and  02  from  the 

first  pair  of  Eqns.  (27)^,  (27^  and  likewise  from  the  second  pair  of 
Eqns.  (27)^  and  (27)^.  Thus  from  Eqns.  (27)  we  can  easily  obtain 
another  set  of  four  equations 


CA^2 

1 

Sh2SA02 

C2P2 

SA62 

1 

Sh2CA°2 

C2P2 

SA*2 

1 

Ch9CA°2 

C2P2 

CA^2 

1 

Ch2SA°2 

C2P2 

C.  ui  * 

hl  1 


Clpl  Shl0;L 


C,_  a 


C2P2  hl 
Clpl 

sh  u 

C2p2  hl 


1 * 

1 ‘ 


(28) 


It  may  be  remarked  that  the  two  sets  of  equations  may  also  be  obtained 

• k 

by  expanding  the  analytic  functions  about  the  stations 

k+1 

and  Y^  and  using  the  continuity  conditions  at  the  three  interfaces 

W , k k+1 

(Y*  + hl^’  ^2  + h2)’  and  ^Y1  + hi')  ’ respectively. 

To  obtain  uncoupled  systems  of  equations,  we  eliminate  and  a 2 

from  the  set  of  Eqns.  (27)  and,  and  from  the  set  of  Eqns.  (28). 

This  elimination  process  thus  leads  us  to  two  pairs  of  uncoupled  equa- 
tions 


(S,  C,  + 0C.  S.  )C.U.  (y,t) 
hj_  h2  hl  h2  A 1 


(ShlSh2+  6Ch1Ch2)SAUl(y’t) 


(ShlCh2+  6Ch1Sh2)CAU2(^t) 


W+  6ShiSh2)SAU2(^t> 


<Cu  Cu  + 6SU  su  )SA0. (y,t)  = 


Clpl  h]  h2  h!  ho  A 1 
1 


0, 

(Cu  Sw  + 0SU  Cu  )CAa1(y,t)  = 0, 

0, 

0, 


Clpl  hl  h2  h!  "2  11 


-(s,  + ec.  a )sAo0(y,t)  = 


(29) 


C2P2  hl  h2  h!  h2  A 2 


(C,  Su  + es,  cu  )cA0,(y,t)  = 


C2P2  hl  h2  h!  "2  A 2 


where  0 = ^2P2^1P1  ’ 
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•k  • 

Let  be  the  determinant  of  the  coefficients  of  and  , 

* . 
and  D2  the  determinant  of  the  coefficients  of  U and  O2  . in  the 

* * 

system  of  homogeneous  Eqns.  (29).  It  is  easy  to  verify  that  = D2  (= 

•k  • • 

say  D)  , and  therefore  each  of  the  analytic  functions  U^,  >°2^=<tl^  ’ 

where  $ is  a collective  symbol,  satisfies  a differential  equation  of 
the  form  6$  = 0.  Explicitly,  the  infinite  order,  linear,  partial 
differential  equation  takes  the  form 


_C2A  - 2<9  + 1/9>  S2h1S2h2  - C2hlC2h2J 


$ = 0 . 


(30) 


By  slight  rearrangement  this  equation  can  be  put  in  the  more  useful  form 


> - * -=J  C2(hl+h2, 


) + if*  - ifc 


/? / 2^hl  h2^J 


4>=  0 


(31) 


We  have  thus  replaced  our  infinite  set  of  governing  equations  (1) 
and  (2)  with  a set  of  four  partial  differential  equations.  From  these 
equations  a series  of  approximate  equations  of  increasing  order  can  be 
obtained  by  retaining  only  a finite  number  of  terms  in  the  symbolic 
operators 


„ - . o,3  " (2A) 2n  , 3 ,2n 

Sa  = cosh  2%  = l j^yr 


'3y 


(32) 


, 2n 


- [2(h  S ±hJ  )] 
C2(h1±h2)  = C°Sh  2(h161±h282)  -37  = I (2n) ! (Tt) 


2n 


For  developing  the  higher  order  equations  we  will  use  Eqn.  (3J),  because 
all  the  symbolic  operators  in  this  equation  are  linear.  However,  we 
will  have  occasion  to  use  Eqn.  (30)  when  we  discuss  the  spectral  function 
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corresponding  to  this  operator  equation. 

The  lowest  order  approximation  to  Eqn.  (31)  can  easily  be  obtained 
if  we  retain  the  first  two  terms  in  the  expansion  of  the  symbolic 
operators  entering  in  this  equation.  Thus,  after  simplification  we 
obtain 


a2 

/hlPl  + h2p2\/hl^nl  + h2^n2^ 

a2" 

U2  ■ 1 

V hl  + h2  /\  hl  + h2  J 

3t2. 

0 < y < 00 , t > 0 


(33) 


This  is  a one-dimensional  wave  equation  and  from  this  equation,  the 
effective  longitudinal  phase  velocity,  density  and  modulus  can  easily 
be  read.  Thus  we  find 

i.  effective  zero-order  density 

(h.^  + h2P 2 ) / (hf  + * 


ii.  effective  zero-order  modulus 

(h1  + h2)/(h1/n1  + h2/n  2)  , 


(34) 


iii.  effective  zero-order  phase  velocity 


(°) 

C2 


(Pj_/p2  + h2/h1)  (n2/n1  + h2/h1)/[C2(l  + h2/h1)]' 


n -1 


(o).(o) 

n/  p • 


The  zero-order  phase  velocity  is  the  effective  celerity  of  the  longitu- 
dinal wave  motion,  with  wavelengths  much  longer  than  the  thickness  of 
the  laminations,  We  may  therefore  call  it  the  long  wavelength  sound 
velocity  of  the  laminated  composite.  This  result  is  well  known  and 
was  obtained  by  Rytov  [6],  and  simultaneously  by  White  and  Angona  [7]. 
The  methods  used  by  these  authors  are  different  from  the  one  developed 
in  this  paper,  and  in  the  second  part  of  this  paper  we  will  show  how 
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the  zero-order  effective  properties  of  a laminated  medium  can  also  be 
obtained  by  using  a multiple  scale  expansion  of  the  displacement  field. 

The  lowest  order  approximation  yields  a non-dispersive  wave  equa- 
tion. The  initial  value  problem  corresponding  to  this  equation  requires 
that  the  Cauchy  data  $(y,o),  $(y,o)  be  prescribed  at  every  point  on 
the  support  0 < y < °°,  t=0. 

Before  developing  higher  order  equations,  we  first  introduce  non- 
dimensional  variables  S and  t , related  to  our  original  variables 
y and  t , respectively,  as  follows 


T : y - A?  , t + (A/C1)t  , 


where  we  have  used  as  the  reference  speed  of  the  lamina  with 

material  properties  n^,p^  . Furthermore,  we  introduce  material,  geomet- 
ric ana  kinematical  properties  in  terms  of  the  ratio  p,  n,  h and  c , 
where 


P = P1/p2  > n = r'1/h2  , fi  = h^/h0  , c = 0^2  . 


In  terms  of  the  new  variables  z,  and  t and  the  parameters  p,  f|,  fi 
and  c , the  linear,  infinite  order,  partial  differential  equation  (31), 
takes  the  non-dimensional  form 


, „ 3 1/  rr  . 1\2  , 2(1  + c/fi)  3 

cosh  2 — - t-  /0  + — cosh  — 

L 35  M /e/  1 + 1/h  8T 

+ cosh  & --  ^ f|  *(5ix) 

4\  /e/  i + i/h 


for  0 < 5 < 00  , t>0  and  where  0 = l/(cp)  = c/fi  . 


Approximate  differential  equations  of  increasingly  higher  order 
are  obtained  from  Eqn.  (37)  by  systematic  use  of  the  expansions  of  the 


4 


symbolic  operators  entering  in  this  equation.  Thus,  partial  differen- 
tial equations  of  fourth  and  sixth  order  are 


i.  fourth  order  approximation 


r'32  2 j|M  2^_  /s4 

9 ~ a 2/  + 4'  \ i 


ii.  sixth  order  approximation 


,4  3 


7/J 


$ = 0 , 


(38), 


L 2 ili , ii  / a4 

2 ~ a 2 + 4 ' 4 

L\3C  3t  / \3C 


4 34  \ 25  /S6  6 36 


4>=  0, 


(38)ii 


where 


(p  + 1/h) (1/n  + 1/h)  ~2 

(1  + 1/h)2 


1 + — + f- 


8 = 


ph 


h / J 


1_  , 2P  1_ 

"2  h -2 


(1  + 1/h) 


(39) 


p + — (l  + n/h)+  r 
h hJ-J 


ri  + i(i  + ±4+  % 

-h  h\  ph/  h J ~2 


(1  + 1/h) 


These  equations  in  their  present  form  contain  fourth  and  sixth  order 
partial  derivatives  with  respect  to  variable  t . However,  we  will 
show  that  we  can  replace  this  set  of  equations  by  an  equivalent  set 
of  equations  which  contain  only  the  second  partial  derivative  with 
respect  to  t , and  thus  requiring  Cauchy  data  pertaining  to  $(y,o) 
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and  4>(y,o)  along  the  support  at  t-0  . 

We  may  now  remark  that  Eqn.  (37)  is  exact  and  replaces  the  infinite 
system  of  Eqns.  (1),(2)  and  interface  continuity  conditions  (3),  (4) 
by  a single  set  of  linear  partial  differential  equations  of  infinite 
order.  If  this  is  exact  then  it  must  yield  the  same  dispersion  spectrum 
as  the  one  we  would  have  obtained  from  our  original  infinite  system  of 
partial  differential  equations  with  piecewise  continuous  and  periodic 
coefficients.  As  a matter  of  fact,  using  Floquet's  theory  of  differen- 
tial equations  with  periodic  coefficients,  Rytov  [6],  and  Lee  and  Yang 
[8],  have  obtained  the  dispersion  equation  for  this  relatively  simple 
problem.  It  will  therefore  be  instructive  to  compare  our  results  with 
the  ones  given  by  these  authors. 

For  dispersive  waves  we  assume  solution  in  the  form  of  sinusoidal 
wave  trains 

$ = i>o  exp  i(q£  - fix)  , (40) 


where  q is  the  Floquet  wave  number,  ft  the  non-dimensional  frequency 
and  the  amplitude  which  is  arbitrary.  We  immediately  notice  that 

if  we  substitute  Eqn.  (40)  in  the  non-dimensional  form  of  Eqn.  (30), 
then  we  get  the  same  spectral  function  as  given  by  Rytov,  and  Lee  and 
Yang.  Having  established  that  our  differential  Eqn. (30)  yields  the 
exact  dispersion  spectrum,  we  now  revert  to  Eqn.  (31)  whose  spectral 
equation  corresponding  to  sinusoidal  wave  trains  is 


cos 


2q  ~ TK8  + 

✓e 


COS 


ILLiiZM  a + kU.  S)\os 
i + i/fi  4\  /el 


2(1  - 8/h) 
1 + l/h 


0. 


(41) 
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As  shown  by  Lee  and  Yang,  the  graph  of  this  function  has  a banded 
structure,  consisting  of  passing  bands  and  stopping  bands  and  the  width 
of  Brillouin  zone  is  tt/2  , Our  aim  now  is  to  show  that  from  Eqn.  (37) 
a system  of  approximate  differential  equations  can  be  derived  whose 
spectral  graph  in  the  first  Brillouin  zone  agrees  with  the  corresponding 
graph  of  the  exact  equation,  up  to  any  desired  degree  of  accuracy. 

Now  according  to  a theorem  by  Gelfand  and  Levitan  [9],  a differ- 
ential equation  can  be  determined  from  its  spectral  function.  If  we 
establish  the  correspondence 

" " 1 k ■ 0 " 1 It  ■ (42) 

then  this  theorem  is  trivially  obvious  in  the  case  of  polynomial  spec- 
tral function.  Using  this  correspondence,  it  is  also  obvious  that  from 
the  spectral  Eqn.  (41)  we  get  the  differential  Eqn.  (37).  Hence,  a 
sequel  to  this  is  the  corollary  that  there  is  a 1:1  mapping  between 
the  (q , ft  )•*•->■  (5 , i)  spaces.  Consequently,  for  those  problems  for  which 
the  dispersion  spectrum  can  be  easily  obtained  by  using  Floquet's 
theory  of  differential  equations  with  periodic  coefficients,  the  corres- 
ponding operator  differential  equation  can  be  immediately  derived  by  < 

using  Gelfand  ~ Levit an  theorem.  Alternatively,  one  may  first  derive 
the  operator  differential  equation  by  using  the  method  of  Lie-series  and 
then  use  this  correspondence  to  determine  the  spectral  function.  Because 
of  the  1:1  mapping  between  the  operator  equation  and  its  spectral  func- 
tion, a power  series  representation  in  the  (q, ft) -space  can  be  translated 
into  a power  series  representation  in  the  (s,i)-space.  Hence,  replac- 
ing the  symbolic  operators  by  their  power  series  representation,  and 
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retaining  terms  of  order  A , gives  us  an  approximate  differential 
equation  of  order  2n.  Furthermore,  the  dispersion  spectrum  corres- 
ponding to  this  system  of  higher  order  equations,  is  just  the  power 
series  representation  of  the  exact  dispersion  equation  up  to  the  same 
order.  Hence,  the  dispersion  equations  for  the  Eqns.  (33)  and  (38) 
can  be  directly  obtained  from  expansions  of  Eqn.  (41) , or  conversely, 
from  the  power  series  representation  of  Eqn.  (41)  , one  can  obtain  a 
system  of  higher  order  differential  Eqns.  (33,38). 

From  Eqn.  (41)  we  readily  see  that  about  the  origin  (o,o) , it  is 

an  even  function  with  respect  to  both  variables  and  has  a functional 

2 2 

representation  of  the  form  F(fi  ,q  ) = 0 . Therefore  by  the  implicit 

2 

function  theorem,  if  3F/3ft  ^ 0 at  the  origin,  then  in  the  neighbor- 

2 2 

hood  of  the  origin  it  has  the  representation  fi  = f(q  ).  Analyticity 
of  the  function  f leads  us  to  the  existence  of  the  power  series  rep- 
resentation 


22  4 4,  6 6 .8  8 

aQq  - eQq  + yQq  - 6Qq  + 


(43) 


Gelfand-Levitan  theorem  now  asserts  that  corresponding  to  this  poly- 
nomial spectral  equation f the  differential  equation  is 


/ 2 a2 

V’l? 


4 34  6 36  8 3 8 

+ 6o  7T  + Yo  76  + 6o  78 

3A  3C  3 e 


••)*  • 


(44) 


If  we  substitute  Eqn.  (44)  in  Eqn.  (37)  , we  get  an  identity.  This 
determines  the  coefficients,  which  are 


] 8 


2 _ i,  2 

a0  - 1/a  , 


4 1 4 4 

eo  ^2  (1  ‘ 6 , 

3a 


6 _ 8 ,,  6 , 6^  2 , .4  , 4,  , . .4,4. 

Yn o (1  ~ Y /a  ) 5-  (6  /a  )(1  - 6 /a  ) , 

45a  9a 


68=-i- 

0 315a 


/■i  r8  . 8.  2 1”.  6,  6.  4,  4.  2 . 4 , 4W1  6,  6. 

j (1  - 6 /a  ) - 45^  /a  )(1  - 8 /a  ) + j(8  /a  )(1  - y /a  ) 


1 4,4  4 4 4 4 

2f  (6  /a  )(1  - 8 /a  )(1  - 58  /a  ) . 


Thus  corresponding  to  Eqn.  (38) we  have  the  equation 


-r  - a 
2 o 


o „ 4 

94 


)<!>  = 0 , 


and  corresponding  to  Eqn.  (38)  ^ we  have  the  equation 


M2 


2 3 2 4 34  6 8 6 

o 2 o A 6 

3£  H 


= 0 . 


As  mentioned  earlier,  we  now  have  a well  posed  Cauchy  problem  requiring 
only  the  two  initial  values  *(^,t  ) and  $(^,tq)  on  the  support 
t = tQ  >_  0 . The  Eqns.  (48)  and  (^8)^,  with  Cauchy  data  given  on  the 
support  t = tq  , determine  4>(4,t)  uniquely  and  therefore  $(c,tq)  and 
$*(4,tq)  can  now  be  determined.  Knowing  3^$(^,to)  for  n = 0,1,2  and  3, 
gives  us  all  the  necessary  initial  conditions,  for  the  initial  value 
problem  governed  by  Eqn.  (38)*  . Analogous  procedure  holds  for  Eqn.  (38). 


This  follows  from  Cauchy-Kowalewski  theorem  as  pointed  out  by  Prof.  R.  P. 


S haw . 
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Corresponding  to  Eqns.  (38/46) ^ the  dispersion  equation  is 


ft2  - o^q2  + B^q4  = 0 , (A7) 

which  leads  us  to  the  phase  velocity  relation 

ft/ q = 1/a  - t~(1  - 64/a4)q2  - ...  . (48) 


This  suggests  that  the  zero-order  phase  velocity  is 


(o) 

C2 


2 2 
q/er 


(49). 


and  the  second-order  correction  to  the  long  wavelength  sound  speed  is 
given  by 


(2)  144  2(0) 
C2  = 3!  (B4/ a4  - l)2  C2 


(49) 


ii 


where  a and  8 are  given  by  Eqn.  (39). 

In  the  case  of  a homogeneous  medium,  n = fi  “ p = 1 and  therefore 
from  Eqn.  (39)  we  find  that  a = B = 1 • Consequently 


(50) 


and  we  get  the  well  known  results  for  longitudinal  v“ves  in  a homogeneous, 
, isotropic,  infinite  one-dimensional  medium. 

Corresponding  to  Eqns.  (38/46)^  the  dispersion  equation  is 

__  
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n2  2 2 4 4 66  n 

A - aq  + g q - y q = 0 


(51) 


which  yields  the  nondimensional  phase  velocity  relation 


ft/q  - l/a.  - ~ (l-64/a4)q“+ 


y(l-Y6/a6)-(e4/a4) (l-64/a4) 


±(1-B4/.4)2 


q4+ 


(52) 


This  suggests  that  the  fourth-order  correction  to  the  long  wavelength 
sound  speed  is  given  by 


(4) 

C2 


1_ 

81 


|(l-Y6/a6)-(64/a4)  (lV/a4)-  |-(l-f34/a4) 2]  " ^ 


(53) 


From  Eqns.  (48)  and  (52)  we  further  see  that  the  limiting  value  of  non- 
dimensional  group  velocity  is 


(dfl/dq)  = 1/a  , (54) 

q-+o 

and  for  long  wavelengths 

d2fi/dq2  « - ^(1  - 64/a4)q  + ...  . (55) 

The  noc-zero  value  of  curvature  for  long  wavelengths  suggests  the  dis- 
persive character  of  the  waves. 

We  may  conclude  this  section  with  the  remark  that  frequency  Eqn.  (41) 
is  of  class  Caj  in  the  neighborhood  of  a point  an^  therefore 

there  exists  a Taylor  polynomial  which  approximates  this  equation  to 
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any  given  desired  accuracy.  From  the  nature  of  the  derivation  of  our 
approximate  differential  equations,  and  as  a consequence  of  Gelfand- 
Levitan  theorem,  it  is  obvious  that  our  approximate  dispersion  spec- 
trum given  by  Eqns.  (47)  or  (52),  are  actually  the  Taylor  polynomial 
approximations  of  Eqn.  (41),  in  the  neighborhood  of  the  point  (o,o). 

By  increasing  the  degree  of  the  Taylor  polynomial,  we  can  obtain  a 
polynomial  approximation  with  arbitrarily  small  error.  This  implies 
that  by  increasing  the  degree  of  approximate  differential  equation, 
we  can  make  the  error  arbitrarily  small.  For  purposes  of  comparison, 
we  have  shown  in  Fig.  2,  over  one  Brillouin  zone,  the  exact  spectrum, 
and  its  various  orders  of  Taylor  polynomial  approximations  as  given  by 
Eqns.  (38)  and  (46).  The  spectrum  marked  and  Lg  , belong  to 

Eqns.  (38)^  and  (38)^,  and  those  marked  and  , belong  to 

Eqns.  (46)^  and  (46) ^ . Spectral  curve  Hg  corresponding  to  Eqn.  (44) 
or  Lg  corresponding  to  the  eighth-order  approximation  of  Eqn.  (37) 
is  also  shown  in  the  figure.  As  mentioned  earlier  in  the  introduction, 
the  spectrum  of  the  eighth-order  theory,  matches  to  a very  high  degree 
of  accuracy,  the  exact  spectrum,  over  the  entire  first  Brillouin  zone. 


Homogenization  and  Multiple  Scale  Expansion 

We  now  derive  the  effective  properties  of  the  laminated,  periodic 
composite  by  homogenization  via  a multiple  scale  expansion.  Before 
doing  so,  we  may  mention  that  homogenization  can  also  be  carried  out 
by  using  functional  analytic  methods  [10,11],  based  on  energy  considera- 
tions and  G-convergence  of  operators  [12] . One  may  also  use  probabilis- 
tic methods  [13],  or  diffusion  approximations  to  random  evolutions  as 
in  [14].  As  a matter  of  fact  these  functional  analytic  and  probabilistic 
methods  provide  a justification  for  the  multiple  scale  expansion  method, 
which  we  now  present. 

Combining  Eqns.  (1)  and  (2)  we  get  the  equation 

T~(n(y)  - p(y)  =0,  0<y<°°,  t > 0 (56) 

3y  \ 3y  / at2 

where  we  now  assume 


n(y)  = n(y  + 2a)  ; p(y)  = p(y  + 2A)  . 


(57) 


The  periodicity  of  the  coefficients  of  the  differential  equation  with 
half  period  A = (h^  + h^) , reflects  the  periodic  nature  of  the  par- 
ticular heterogeneous  medium  under  consideration,  and  we  wish  to  study 
Eqn.  (56)  as  A -+•  0 . 
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and  rewrite  Eqn.  (56)  as 


(l(x)  UL(x,y;t)j  - p(x)  UL(x,y;t)  . 

3 1 * 


In  this  equation  the  superscript  e indicates  the  dependence  of  U on 
e . We  want  to  study  the  behavior  of  Eqn.  (59)  as  e ■*  0.  Keeping  y 
fixed  and  letting  e approach  zero,  we  find  that  n(x)  and  p(x) 
become  more  and  more  oscillatory.  Therefore,  the  study  of  Eqn.  (59) 
as  e ->  0 , does  provide  us  with  information  about  Eqn.  (56)  as  A -*■  0. 

In  view  of  Eqn.  (58) 


3/3y  = 3/3y  + ~ 3/3x  , 


and  therefore  Eqn.  (59)  can  be  rewritten  as 


2t  j e i p 

+ - n(x)  + — n(x)  — — + — — n(x)  - — = p(x)  

3y2  t L 3y3x  3x  3y  J e2  3x  9t2 


32UE  , 3 


3U  1 . 1 3 


We  now  assume  that  U (x,y;t)  has  a power  series  representation 


UE(x,y;t)  = l en  U (x,y;t)  , 


where  it  is  assumed  that  each  of  the  functions  U , n = 0,1,2,...  are 

(n)  (n) 

2A-periodic  with  respect  to  x ; i.e.,  U (x  + 2A,y;t)  = U (x,y;t). 
Substituting  this  series  representation  in  Eqn.  (61)  and  equating  like 
powers  of  e , we  get  the  system  of  equations 


f 


Equation  (63)  ^ for  n = 0 , has  a periodic  solution  with  respect 


to  x iff 


^ Jo)  Jo) 

(.)  ^ ^ -nC) 

L 3t2  3y2 


XD  (D- 

92U  3 n 3D  , „ . 
— — — - — n(x)  — dx  = 0 ; 
3y3x  3x  3y 


since  the  remaining  term  when  integrated  over  one  period  vanishes. 
Substituting  Eqns.  (64)  and  (66)  in  (68)  and  simplifying,  we  get 


r \ *o)  (o)  i°> 

t ,wa  4 df,«d,  a . £S  . „ , 


where 


x+2A 

(°)_  1 f 

p E 2A  p(x)dx 


is  the  average  density  of  the  material  over  one  period.  To  find  the 
value  of  the  integral  in  Eqn.  (69),  we  make  use  of  Eqn.  (67).  Integrat- 
ing this  equation  once  we  get 


n(x) (1  + f^(x))  = n , 


where  ri  is  the  constant  of  integration.  To  determine  this  constant 
we  integrate  Eqn.  (71)  over  one  period,  and  using  the  periodic  property 
of  function  f^(x)  , we  get 


T^T  = Ia  j (l/n(x»dx  , 


which  is  the  average  of  q (x)  over  one  period.  Now  substituting 


Eqn.  (71)  in  Eqn.  (69),  we  get  the  "homogenized  equation" 


26 


(o) 


= 0 , 


(73) 


3y^ 


3t" 


which  is  the  zero-order  approximation  to  Eqn.  (56)  or  (61).  This  is 

a one-dimensional  wave  equation  in  a homogeneous,  isotropic  medium 

(o) 

whose  effective  zero-order  density  and  stiffness  are  given  by  p 

(o)  (o) 

and  q , respectively.  In  the  case  of  a bilaminate  composite,  p 

(o) 

and  q calculated  from  Eqns.  (70)  and  (72)  agree  with  the  results 

given  by  Eqns.  (34).  respectively.  This  verification  of  the 

i,n 

results  obtained  earlier  by  using  Lie-series,  provides  us  with  another 

strong  justification  for  the  multiple  scale  expansion  just  carried  out. 

(2)  / „ 

To  determine  U we  again  consider  Eqn.  (63) ^ for  n=0.  Making  use 
of  Eqns.  (64), (66)  and  (73)  we  get 


3_ 

3x 


(2)  (1) 

, . 3 U 3 U 

n(x)  + ^ 


(1),  . 32U  _ 

p (x)  — 2 0 

3t 


(74) 


where 


(1)  o 

p (x)  - p(x)  - p , 


(75) 


which  is  the  deviation  in  density  from  the  effective  density  p . 

We  assume  that  the  solution  of  Eqn.  (74)  , which  is  2A-periodic  in 
x is  of  cne  form 


«) 

f2(x)  + U 

3y 


2(o)  (o) 

3 U , . n3U 


(76) 


= f2(x) 


o 

+ U . 
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Now  substituting  Eqn.  (76)  in  (74)  and  making  use  of  Eqns.  (71)  and 

(2) 


(73),  we  find  that  for  U to  be  the  2A-periodic  solution,  f,,(x) 
must  satisfy  the  equation 


j (o)  o 

^[n(x)(f’(x)  + fx(x))]  = (P(x)  - p )Cq  , 


(77) 


2 (o) (o) 

where  Cq  = n/  p . Integrating  it  once  we  get 


1 2 (1)  (1) 
f2<*>  ’ iW(C.  5<x)  + " > - 'l«  • 


(78) 


(1) 


where  n is  a new  constant  of  integration.  Integrating  this  equa- 
tion over  one  period  and  making  use  of  Eqn.  (72)  , we  get 


„2  x+2A 

c r (1) 


o 

2A 


(l)(o)  1 

p (x)dx  + n/  n 


x+2A 


2A 


f1(x)dx 


(79) 


However,  from  Eqn.  (71) 


f1(x) 


(o) 

n 


dx 


(1) 

o n(x) 


(80) 


where 


(1)  (o) 

1/  i)(x)  = (l/n(x)  — 1/  n ) , 


(81) 


which  is  the  deviation  in  compliance  from  the  effective  compliance. 


Now  substituting  Eqn.  (80)  in  (79),  we  find  that  the  effective  first 


(1) 


order  stiffness  n is  given  by 
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n/  n = IT  I j 


(o)  / x 


dx 


d£ 


(1) 

o n (c) 


x+2A 

r (1) 


j 


(o) 

( 2(x)/  p )dx 


(82) 


This  process  can  now  be  continued  and  effective  properties  of 
higher  orders  can  similarly  be  determined  from  Eqns.  (63)  , n » 2,3,.. 
We  now  note  that  when 


(1)  3U}  ° 

u = fi(x)i!r  + u * 


2 (°)  (°) 

(U)-  £,<«)  i-S-  ♦ £,(,)  -f-  * 8 , 

2 3y2  1 3y 


(83) 


(n)  f r , 3n(S}  . . ( . 3n'l(?  . 

U " fn(x)  — f <x)  —^+  * 
3y  n-1  3y 


<°> 

• + fl(x)  13T+U  > 


the  power  series  representation  for  displacement  field  takes  the  form 


U£(x,y;t)  = (1  + ef1  + e2f2  ~ + •••)V(y,t)  , 

y 3y 


(84) 


where 


V(y,t)  = 8(1  + c + e2  + ...)  - U 


(1-e)  ’ 


(85) 


and  satisfies  the  "homogenized  equation"  (73). 

e ° 

It  can  be  shown  (cf.  [14])  that  U converges  weakly  to  U in 


L2((P,«>),  H^(O.oo))  and  that  UE(x,y;t)  - (1  + ef-^x)  |^-)V(y,t) 

converges  strongly  to  0 in  the  same  space,  (cf.  [10]).  Here 

H1(0,“)  = {f (y)  : f and  e L2(0,“),  f(o)  = f (“)  = 0}  . 
o oy 
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General  Remarks 

The  general  solution  to  Eqns.  (1)  and  (2)  is  a quasiperiodic  Floquet 
wave  of  the  form 

U(x,t)  = v (x)  exp [ i (qx  - ftt) ] , (86) 

where  the  amplitude  function  v is  2A-periodic  in  x,  i.e.,  v(x  + 2A)  = 
v(x) . 

As  was  shown  in  part  one,  the  Lie-series  approximation  amounts 
to  assuming  that  the  frequency  ft  can  be  expanded  in  the  power  series 
(43).  Thus,  the  Lie-series  approximation  is  a long  wavelength  approxi- 
mation. 

The  approximation  obtained  by  homogenization  in  part  two  corres- 
ponds to  an  expansion  of  the  amplitude  function  as 

v (x)  = 1 + iq  v1  + (iq)2v2  + ...  . (87) 

Note  however  that  the  wave  number  q in  this  expansion  corresponds  to 
that  of  a wave  propagating  in  a homogeneous  medium,  the  properties  of 
which  are  the  effective  nroperties.  This  approximation,  in  which  q 
is  of  the  same  order  as  ft  but  differs  in  order  from  the  q entering 
in  the  Lie-series  approximation,  is  therefore  a low  frequency  approxi- 
mation. 

The  approximation  suggested  by  Kohn  in  [4]  is  one  in  which  both 
frequency  ft  and  amplitude  v are  expanded  in  power  series  in  wave 
number  q . It  is  thus  a long  wavelength-low  frequency  approximation 
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as  pointed  out  by  Kohn  himself. 

The  merits  of  each  approximation,  which  will  dictate  the  choice 
of  one  approximation  over  the  others,  depend  upon  the  region  of  interest 
in  the  (fi,q)  space  as  well  as  the  nature  of  the  problem  under  investi- 
gation (e.g.,  transient  versus  steady-state  problems). 

The  way  in  which  microstructure  is  taken  into  account  in  the 
homogenization  expansion  is  worth  a few  comments.  Indeed,  f^  intro- 
duces a correction  for  variable  stiffness  only,  (cf.  Eqn.  (79));  a 
correction  for  variable  density  first  occurs  through  f^  , (cf.  Eqn.  (78)). 
Both  stiffness  and  density  enter  into  the  higher  order  correction  terms 
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